Spin Glass and ferromagnetism in disordered 
Cerium compounds 

S. G. Magalhaes*'^ F. M. Zimmer^ P. R. Krebs^ B. Coqblin^ 

"'Lab. de Mecdnica Estatistica e Teoria da Materia Condensada, 
Dep. Fisica, UFSM, 97105-900 Santa Maria, RS, Brazil 
^ Instituto de Fisica e Matemdtica, Universidade Federal de Pelotas, 
Caixa Postal 354, 96010-900 Pelotas, RS, Brazil 
^ Laboratoire de Physique des Solides, Universite Paris-Sud, 
batiment 510, 91405 Orsay, France 



Abstract 

The competition between spin glass, ferromagnetism and Kondo effect 
is analysed here in a Kondo lattice model with an inter-site random 
coupling Jij between the localized magnetic moments given by a general- 
ization of the Mattis model |15| which represents an interpolation between 
ferromagnetism and a highly disordered spin glass. Functional integral 
techniques with Grassmann fields have been used to obtain the partition 
function. The static approximation and the replica symmetric ansatz have 
also been used. The solution of the problem is presented as a phase dia- 
gram giving T/J versus Jk/J where T is the temperature, Jk and J are 
the strengths of the intrasite Kondo and the intersite random couplings, 
respectively. If Jk/J is small, when temperature is decreased, there is a 
second order transition from a paramagnetic to a spin glass phase. For 
lower T/J, a first order transition appears between the spin glass phase 
and a region where there are Mattis states which are thermodynamically 
equivalent to the ferromagnetism. For very low T/J, the Mattis states 
become stable. On the other hand, it is found as solution a Kondo state 
for large Jk/J values. These results can improve the theoretical de- 
scription of the well known experimental phase diagram of CeNii-^Cux 
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1 Introduction 

The properties of many cerium or uranium compounds are well described by 
the Kondo-lattice model, with strong competition between the Kondo effect on 
each site and the Ruderman-Kittel-Yosida-Kasuya (RKKY) interaction between 
magnetic atoms at different sites. The role of disorder has been studied in dis- 
ordered alloys containing cerium or uranium and different theories have been 
proposed. In the Kondo disordered model (KDM) 0|2|, disorder produces a 
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broad distribution of Kendo temperatures and can be responsible for the de- 
viation from the Fermi Liquid behavior found in some heavy fermion systems. 
Another theoretical approach is the Magnetic Griffths phase i3j, where fluctu- 
ations of the magnetic clusters can produce GrifRths-McCoy singularities close 
to a Quantum Critical Point (QCP). On the other hand, we have studied within 
a mean field approximation the phase diagrams observed in disordered heavy 
fermion systems showing Kondo, spin glass and magnetically ordered phases 
01 13 m and we will discuss these models later on. Earlier studies have also 
suggested that a spin glass transition near a QCP could lead to a Non Fermi 
Liquid (NFL) behavior [Jj. Our paper is an attempt to improve the theoret- 
ical description of the spin glass-Kondo-ferromagnetic competition in order to 
obtain a better agreement with the experimental situation of disordered cerium 
or uranium heavy fermion systems. 

Spin glass and Kondo state have been observed together in several Cerium 
alloys hke CeNii-xGux IHUHl, Ce2Aux-xCoxSi3, ^21 and in some disordered 
Uranium alloys such as UCu^-xPdx 53 or Ui-xLaxPd^Al^ jj^. The first case 
has been studied by bulk methods (see Refs. [5] 1^1) and more recently by /iSR 
spectroscopy ^U] which gives local information about the spins configurations. 
The bulk probes have shown a presence of antiferromagnetic phase for low Ni 
content (for instance x = 0.9). In the region x < 0.2, the Kondo effect becomes 
important producing magnetic moment reduction. For 0.8 ^ x > 0.4, the same 
bulk probes have shown a presence of the spin glass like state intermediate in 
temperature between a ferromagnetic order (at lower temperature) and param- 
agnetism (at higher temperature). However, the /xSR spectroscopy has shown 
in the region 0.8>x>0.4a scenario favoring the presence of a inhomogeneous 
"cluster spin glass" (or called equivalently "cluster glass" in Ref. PO]) rather 
than a standard spin glass. Quite recent measurements on the specific heat 
[TT] has confirmed the emergence of a spin glass-like state and a percolative 
evolution to a ferromagnetic order at low temperatures. 

There has been a theoretical attempt to build up a global phase diagram 
based on a Kondo lattice model with a random Gaussian inter-site coupling 

among the localized spins with mean 2Jq/N and standard deviation \J 8J^/N 
{N is the number of sites). The spins operators have been given as bilinear 
combination of creation and destruction fermionic operators. The partition 
function has been found using path integral formalism within the static ap- 
proximation and replica symmetry ansatz The results have shown that 
ferromagnetism, spin glass and a mixed phase (a solution with non-zero magne- 
tization below the Almeida-Thouless line) have been obtained for small Jk/J 
values and Jq/J > 1.46 while a Kondo phase is obtained for large Jk/J values. 

However, the calculated spin glass freezing temperature (Tf) is lower than 
the Curie temperature (Tc) in this highly frustrated model. Even the transition 
temperature to the mixed phase is always below the onset of the ferromagnetic 
order at Tc- Thus, our previous model 5^ gives a ferromagnetic transition 
temperature (Tc) above the spin glass transition temperature (T/), in contrast 
with the experimental situation of CeNii-xCux jSl 1^ IIOL alloys where the 
ferromagnetic phase is always the lowest one. That would be a clear indication 
that Gaussian distributed random couplings as in the Sherrington-Kirkpatrick 
(SK) model is not adequate to describe the frustration present in that alloy. 

One important point in the set of experimental works has been to clarify 
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the effect of the disorder in the CeNii-^Cux [3 El El E]- When Cu is ran- 
domly replaced by Ni in that alloy, not only the cell volume is modified, but 
also the number of conduction electrons which makes the competition between 
the RKKY and Kondo effect complex for that particular alloy jjl]. As a con- 
sequence, one localized spin at any site can be subject to a set of effective local 
magnetic fields, which results in the complicated combination of states previ- 
ously cited. 

In the Mattis model J 5^, which has been proposed as a solvable model to 
the spin glass problem, the bonds joining the localized spins have been defined 
as separable random variables ^i. This model could allow to gain some insight 
about the local effects of disorder as long as it would be possible to construct 
local applied fields dependent on the random variable ^i. Unfortunately, at 
zero magnetic field, a gauge transformation of the Ising spins classical vari- 
ables leads Mattis model free energy to behave rather as the usual ferromagnet 

I17j . Therefore, this model is trivially disordered in the sense that it is 
unable to produce the essential component of the spin glass which is frustration. 
Nevertheless, the generalization of the Mattis model TF, has proved to be 
an interesting alternative. In this model the coupling between spins are given 

by 

where = ±1 (/x = l,2,...,p; i = l,2,...,iV) are independent random dis- 
tributed variables. For the classical Ising model, ii ^ ~ v — 1, the original 
Mattis model ^HI is recovered. However, if J^j, = JS^j^^, and p=N with the 
random variables having mean zero and variance one, in the limit of N large, 
Jij tends to a Gaussian variable with mean zero and variance as m 

the SK model ^5]. Therefore, we can consider this model as an interpolation 
between ferromagnetism and highly disordered spin glass 122' . 

An important particularity of this model (see Eq. Ql) is J^i, — JS^^i, which 
has been used in a different context, i. e. the statistical mechanics theory of 
complex systems |22) using classical Ising spins. In this problem, randomness 
effects can be better understood at T = temperature when the local field 
applied in a particular spin given by hi = X^i^^j Jij^j is analysed |21| . In the 
state Si = (choosing J = 1), the local field becomes hi + 5i), where 

5i is a random variable with variance < <5| >{= Two situations can be 

identified when N ^ oo {N is the number of sites). If p is finite, the spin is 
perfectly aligned with ^f. However, if p increases linearly with N {p = aiV), 
the term Si can become important and the alignment can be destroyed. This 
random component of the local field can be a source of frustration and, in that 
sense, the ratio yN/p is the analog of Jq/ J |221- When temperature is turned 
over, there is an additional mechanism to avoid the aligning. 

Actually, the thermodynamics of the generalized classical Mattis model can 
be described by a mean field theory I23L I24L 171] in terms of the parameter 
a = p/N . For a particular value called Oc, two clearly distinct regimes can be 
identified. When a > Oc, the frustration is dominant below a certain tem- 
perature Tf leading the problem to a spin glass behavior. Nevertheless, when 
< a < Cc, a much more complex scenario can appear depending on the tem- 
perature. For instance, it can be found as solution not only a spin glass phase, 
but also Mattis states (which corresponds to the stable aligning between ^'s and 
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spins) with a first order transition between them. These Mattis states have the 
same thermodynamics as the ferromagnetic phase ^2 "ST, "ST, '7^ . 

The mean fieid description of the previous situation introduces the order pa- 
rameter — J2i < Si > which gives a measure of the difference between 
the configurations of the set {^f } and the spins. In that approach, one partic- 
ular solution can be, for example, only with a possible non zero value, while 
the remaining m'^ {fi ^ 1) are of order 1/\/N. This choice corresponds to the 
situation where the spins can be perfectly aligned only with S,}. Nevertheless, 
these m'^'s (/x 7^ 1) still have a role in the problem yelding a possible spin glass 
solution in the problem depending on the temperature and, particularly, on the 
parameter a. Therefore, this approach would be mathematically convenient to 
apply in a problem where we would be interested in getting control about the 
degree of frustration. 

In this work, we consider the Kondo lattice model with a random intersite 
interaction between the localized spins where the coupling is given by Eq. 
(^. That would allow us to investigate the competition between Kondo effect 
and magnetism combining the approach of Ref. 0] with Ref. Therefore, 
the degree of frustration a = p/N is a new parameter which together with Jk / J 
[Jk is the strength of the intra-site Kondo coupling) constitute the parameter 
space where the solutions can be located. In the Ref. 0], the solutions for the 
order parameters have been found only in the limit of strong frustration (the 
random Gaussian coupling) which corresponds to the a » ac situation. We 
will show a more complex situation in the limit of weak frustration. For small 
Jk / J , there is an intermediate spin glass phase between paramagnetism and the 
region where there are Mattis states which corresponds to the ferromagnetism. 
Furthermore, the transition from the spin glass to the ferromagnetism is first 
order. Therefore, there is a large region in temperature where the ferromagnetic 
solution is thermodynamically metastable. It becomes stable at very low 
temperature. For large Jk/J, we get a Kondo state, as already discussed in 
Ref. II. 

It is not obvious that the properties of the classical Mattis model and its 
generalization can be extended to the quantum version of these models. How- 
ever, it has been shown that the long range quantum Mattis model has the same 
qualitative behavior of his classical counterpart |^. Moreover, in the present 
fermionic spin glass mean field approach with the static approximation, which 
is reliable at high temperature, the spin variables have the essential features of 
classical ones 12Z| ■ 

It should be remarked that the extension of the method given in Ref. 
to the present fermionic problem is not straightforward. In fact, the model 
developed here is different from that one introduced in Refs. Indeed, 
both models study the competition between Kondo effect and spin glass, with 
eventually an additional magnetic phase. But, the description of the spin glass 
onset is different in the two cases. In the present one, the parameter a can tune 
the spin glass or the ferromagnetic component of the long range internal field 
which has a completely different dependence on the spin glass order parameter 
as compared with Refs. Actually, for large values of the parameter a, 

the Gaussian distribution of spin couplings is recovered. That corresponds to the 
strong frustration situation studied in Ref. j4,. Thus, the theoretical description 
is not simple, but the present model is able to give a more local description of the 
problem. Finally, as it will be discussed later on, we will obtain a ferromagnetic 
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phase below the spin glass solution in terms of temperature, in contrast to the 
results of the previous model and in improved agreement with the experimental 
phase diagram of CeNii-^Cux alloys. 

On the other hand, we can say that the experimental situation of CeNii^xCu, 
or similar alloys are also very complex. In fact, the so-called spin glass phase is 
more exactly a "cluster spin glass" which tends to a real ferromagnetic order 
by a percolative process when temperature decreases. 

The outline of the paper is the following. In section II, the model is in- 
troduced and developed in order to get the free energy with the relevant order 
parameters for the problem. The results obtained are presented and discussed 
in section III, and we finished it up with a conclusion in the last section. 



2 General Formulation 

The model is the Kondo lattice used previously to study the competition be- 
tween spin glass and Kondo effect W . The Hamiltonian is given by 

H = Y, ^knl, + E ^0"- + + Jk J2 ISpi + Sjjt] (2) 

where the sum is over the N sites of a lattice. 

The term Hsg corresponds to the intersite interaction between localized 
spins, thus 

Hsg = E ■^v^h^fj- 

The random coupling Jy present in the previous equation is given by Eq. 
^ with Jfj^i, = ^5^1, 122], where — ±1 are random independent variables 
which follow the distribution: 

i^(CP = l/2%.+i + l/2^j;.-i. (4) 

The spin operators in Eq. Q are defined (see Refs. O ^ E|) as bihnear 
combinations of the creation and destruction operators for localized (conduc- 
tion) fermions f}^, f^i (dj^ dii) with the spin projection a =t or J,: S^^ = A^/ii; 
SJ, = 4 /it; s+ = dl^d,i; s- = dl^d,y, 

The chemical potential for the localized and conduction bands are ^/ and 
fic, respectively. As it has been done in Ref. [HI El El, the energy is referred 
to jjbf and e/c is referred to /ic. 

In the functional integral formalism, the partition function is expressed using 
anticommuting Grassmann variables (picr{T) (related to conduction electrons) 
and il^icr {t) (related to the localized electrons) j4j as 

Z = ( D{r^)D{v*^) exp [AsG + Ak +Ao(V'*, i>) + A^^^*, ^)] (6) 
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where in the static approximation ^ 

ij.a.u) 

ia ui jcr ^' 

Astat _ \ " T O O. 
^SG — / /JijJiOj- 



A stat 



(7) 
(8) 

(9) 
(10) 



with 



= ^EE'^'^«o('^)V'*o(w)- 



(11) 



a=± 



In Eq. (|11() . the Matsubara's frequencies are given, as usual, by = (2m + \)-k 
with TO 0,±1,±2, .... 

The problem is treated closely to the mean field approximation of Ref. 
Therefore, the Kondo order parameter and its conjugate are introduced 
using the integral representation of the delta function as 



exp < i E 



(12) 



where the presence of the order parameter A(j is related to the formation of 
d — / singlet throughout the whole lattice. This mean field order parameter is 
presently recognized to provide good description of the Kondo effect on each 
site EHl 

Therefore, the resulting partition function becomes: 



Z = y JldA^dA^ exp 



(13) 



where 



+PJk E [ E ^,l('^)^..(^) + A. E ^U^)VjA^) ] } ■ (14) 



In fact, in this work the Kondo order parameter is considered Ao- = A 

The integration over the Grassmann fields ip* and ip in Eq. H14|l can be 
performed, which results in: 



Zstat 

~Z^ 



(15) 



6 



where 



with 



The Fourier transform of Green's function 7^ {lj) in Eq. (|17|l is 

1 



(16) 
(17) 

(18) 



In order to introduce the proper set of order parameters in our problem |23|. 
the action Af^ in Eq. I|l()|l is written using Eq. Q to give 



A 



SG 



f-L—1 i 



(19) 



where Si has been defined in Eq. Hll|l . 

The free energy can be obtained following the replica method, 

/?/ = 2(3Ju\*X - lim ^ mi\at))i - 1) , 

n— »0 iV n 



(20) 



where ((• • • ))c the averaged over ^'s. 

The fundamental issue consists in the evaluation of the quadratic form 
present in the first term of A^^q . First, the sum over /i is separated in two 

parts El: = ^=, + n;l • 

It is possible to linearize the problem introducing n x p auxiliary fields 

and mj^ (a is a replica index) which correspond to the parameter discussed in 

the previous section. Therefore, 



exp (A; 



stat \ 
SG ) 



exp 



N 



/oo ( «-l r -, 1 1 "I poo 

^Dmjjexp /3j7v^^^--(m:;)2+ -^er^fK y_^c 



m' X 



(21) 



where Dm^'^''^ = [] JJdm^^'V V2^- 

In this work, the structure of solutions for auxiliary fields rn^'s and mj^'s is 
the same as in Ref. We assume that the relevant contributions come from 
m'^ which are order unity while m'^ is of order l/v^- Therefore, the average 
over the p~ s independent random variables can be done using Eq. which 
results in: 



((exp 



/3'^EE(E^r^° 



))? = exp 



^£ln(cosh(/3J^5rmJJ)) 



(22) 
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The argument of the exponential in the right hand side of Eq. 122|l can be 
expanded up to second order in The result is a quadratic term of the spins 
variables iSf in the last exponential of Eq. 1)21(1 . This term can be linearized by 
introducing the spin glass order parameter qap using the integral representation 
of the delta function as we have done with the Kondo order parameter, so 



/OO 1 



(23) 



where Qafs is equivalent to the usual spin glass order parameter introduced in 
the classical SK model which gives the transition to spin glass phase when 
= 0. 

Therefore, after rescaling — ^ m'^/^/N, the exponential involving in 
Eq. (|21|) can be written as: 



exp • 



f ^ r 1 1 



lj,— s a 



1 \ (3 J ^ 

Y 1 '^^^pdfap) exp <^ — — ^ m^Aa^mj; 



af3 



where the matrix element 



Aq/3 — Safi — PJqa/S 



(24) 



(25) 



Introducing Eqs. ((24)1 and ((25(1 into Eq. II21(I . the fields can be integrated 
to give: 



((exp^|*g*))e=exp 



2N 



N 



Q=li = l 



Dm":. 



s-l 



exp<]/3JAr^^ -^(m^)2+ l(^^r^f)< 



+ 00 



(IT '^'^"^^"^ ) 



af3 



2tt 



1 ^ 1 

a/3 i=l 



Therefore, the averaged partition function (see Eq. 1120(1 ') is obtained from Equa- 
tions 1(15117(1 and ^ as 



+ C)C 



-dqapTap 



2n 



exp{-/3iV[^5](m:;)2 + I^TrlnA- 



(27) 
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where 

n = f Diri') exp{-z ^ r^p{^ ^ SfS^) - 

E E E v-*. j^k^'(^)V',.a(^)} (28) 

ij (J a u 

The free energy (see Eq. EOl is evaluated at the saddle point and is given by 
the condition that the first derivatives of integrating variables are zero. There- 
fore, for instance, we have 

^2j2 /32j2 

and 



a/3 



^(to^to^)) = — — p a 7^/3. (30) 



The parameter has been found as in the section I. The problem is treated 
assuming the replica symmetric ansatz, therefore the order parameters are qaf3 = 
<?, Qaa = q, TaoL = Tap = r and = to"^. On the other hand, the trace of the 
matrix A is obtained in terms of its eigenvalues (see Eq. I25|) . In consequence, 
the free energy can be written as 

+ - \n[l - (3J{q- q)\ H — r q — r q - 

lim47M(M^^"^MA|)))d (31) 

where a — p/N. 

The sum over replica index produces quadratic forms into the function 
r2(r, r, m'^ , | A |). This term can be linearized by a Hubbard-Stratonovich trans- 
formation ^ 121 where new auxiliary fields are introduced in the problem. 
Therefore, we have: 

f](r,r,mM A|) = 

+00 N „ + oo N n „ 

n^^w n n ^< / ^(v^*'^) 

^^pE E E ^4ac.{^)Gij\L0)ij,.o.{L0)] (32) 
ij a.cr uj 

with Dz = -^e-^'/2 and 

G^'(^) = 9-/{^) - 'T(7^r(r,r) + pJ^^-mnS., (33) 
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where g.j^j^{uj) is given by Eq. lfTC|) and the local spin glass component of the 
internal field 



h, (r,r) = y/(3Ja[(3J{r - r) - 1]< + f3Jy/^z,. (34) 
The functional integral in Eq. (p?^ can be performed (SlEIEli so we get 

N 



/ + 00 ^' 
(11^ 

Af n 



^^ln[dctG,r.i(c.)] 



(35) 



The local field (r^r) applied in the n replicated lattices is, in fact, asso- 
ciated with the replica diagonal and non-diagonal spin glass order parameters. 
The presence of the local field is the fundamental technical issues which must be 
solved in order to proceed the calculations. Therefore, we consider the following 
decoupling 

ln[detG-.i(c.)] « l^ln[detr,v.(c^,7^r(^,^),^r)] (36) 

i 

which means that a constant field is applied in a fictitious Kondo lattice, so 
the problem can be solvable by a Fourier transform. The particular form of 
the decoupling is also useful because it allows the use of self-averaging property 
W^ifi^i) = ((/(C)))c which is valid in thermodynamic limit for finite 
s — 1, the upper value of v. Actually, from now on, it is assumed that s = 2. 
The self-averaging in Eq. H36|l allows us to drop the site index i in Eq. . 

The resulting sum over k in Fourier transformed Tk^ai^, h (r, r),^) can be 
replaced by an integral using a constant density of states for the conduction 
electrons p(e) ~ for —D < e < D. On the other hand, the sum over the 
Matsubara's frequencies in Eq. (|35|l can be performed following closely the 
procedure given in Refs. Finally, the free energy can be found as 

2 72„ a2 j2„ /•+00 p+co -, /•+/3D 



Dz{{ln{ / i?u;exp[— / dx 

-00 J-oc J-f3D 



2 ^9 2 — ^ 1 ^ I ^•^vU'^i / ^''^^^Pl'^T^ I dx X 

ln[2 cosh( '^ + )] + 2 cosh( VA(r, r, 0^ + \ A l)^]}))^ (37) 



where 



and 



Air,T,0-'-^^ (38) 



h{r, r, = ^Jl3Ja[l3J{r - r) - l]w + ^J^/^z + (3Jm^^ (39) 



is a long range internal field composed by two parts, a spin glass one already 
introduced in Eq. (|34|l and other one associated with the order parameter 
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Figure 1: The phase diagram with T/J versus Jk/J for a = 0.02 showing 
the phases SG (spin glass), FERRO (ferromagnetism) and KONDO (the Kondo 
state). In the region SG+Ferro, there is coexistence between SG and FERRO. 
The dotted Hne means the "pure" Kondo temperature. 

m^. This result can be compared with Ref. [S]. In that work, a random 
Gaussian inter-site coupling among the localized spins has been used. The 
equivalent internal field found there can be also decomposed in two parts: 
a ferromagnetic one, and a spin glass term associated with q and the static 
susceptibility x — P{q~ <?)• However, the dependence of the internal field with 
q and q is entirely distinct in the present work as it will be shown below. 

The coupled set of equations for the order parameters can be found from 
Equation H37() using the saddle point conditions. In particular, there is a rela- 
tionship between the following order parameters: 

[l-f3J{q-qW ' 



pj[l-(iJ{q-q)] • ^ ' 

If we replace Eqs. H40I41II in Eas. H37l39(l . the theory becomes explicitly de- 
pendent on the q and q. Therefore, the minimum set of order parameters to 
be solved in order to obtain a global phase diagram is the spin glass order pa- 
rameters q, q (related to diagonal matrix qaa), the Kondo order parameter | A | 
and the = m. The average over ^ in Eq. (|37|l can be now performed using 
the parity properties of the functions dependent on z and w. Therefore, the 
dependence on ^ can be dropped. The remaining set of coupled equations for 
the order parameters are given by the corresponding saddle point equations. 



3 Discussions 



In this paper, it has been investigated the competition between ferromagnetism 
and spin glass in a Kondo lattice model with a random coupling between 
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Figure 2: The behavior of the order parameters q, q and m showing the second 
order and the first order transition for Jk ~ 5J. 

localized spins. The coupling is constructed (see Eq. |^) as a product of 
two random independent variables where i = 1...N (N number of sites) and 
fj, — l...p which is a generalization of the Mattis model |15| for random magnetic 
systems. In the strongly frustrated limit p — > cx3 (when N oo and p/N — 1), 
it is recovered a spin glass solution as S-K model The problem is solved 

using functional integral formalism, the static aproximation and the replica 
ansatz. The order parameters are obtained by combining methods proposed in 
Refs. [22] which allow us to introduce an aditional parameter a = p/N to 
control the level of frustration. 

The numerical solutions of the order parameters q, q, m and | A | (from 
now on I A |== A ) give as solutions the following thermodynamic phases: (a) a 
spin glass where q with A = and m = 0; (b) a Kondo state with A 7^ 
and q = and m = ; (c) ferromagnetism which is given by the existence of 
the Mattis states described by m 7^ and q ^ with A = 0. The solutions 
are displayed in diagrams T/J versus Jk/J for several values of a, where T 
is the temperature, Jk and J are the strength of the intrasite Kondo coupling 
and the intersite random coupling between localized spins (see Eqs. ©-(Ell), 
respectively. Therefore, for a given Jk/J, it is possible to probe solutions for 
the order parameters equations q ,q, m and A in several situations ranging from 
weak frustration to strong one just by varying the parameter a. 

The physical origin of the phases discussed in the previous paragraph can be 
understood from the model introduced in Eqs. (2)-(3) in which there are two 
interactions. The first one is on site Kondo coupling while the second one is the 
disordered coupling between localized spins. Therefore, it is possible to identify 
several energy scales in the problem as T, Jk and J. When temperature is 
high enough, there is only paramagnetism. As long as the Kondo energy scale 
becomes dominant in relation to the remaining ones, the emerging ordering is the 
complete screening of the localized spins in the whole lattice due to the Kondo 
effect. However, in a certain range of Jk/J, there are two possible magnetic 
orderings. In the spin glass one, the competing interactions between the spins 
can give rise to frustration in which there is a large number of degenerate states 
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Figure 3: The free energy as a function of temperature for a — 0.02, a = 0.04, 
and a = 0.15 showing the region (between Ti and T2) where there are multiple 
solutions for the order parameters. 

for the spins configurations. Therefore, there is no long range order correlation 
among spins orientations due to the presence of frustration which leads the 
spins to be frozen in random orientations. In contrast, in the ferromagnetic 
regime, the Mattis states correspond to the situation in which the spins stabilize 
aligned with the ^^'s due to the low frustration level. The prevalence of one or 
other regime discussed above depends on the parameter a which controls the 
frustration level for the random coupling given in Eq. (1) as well as T. 

It should be remarked that previous investigation, using the same model, has 
adopted the standard Gaussian random coupling ^ (as the S-K model) for 
the intersite interaction between localized spins. Nevertheless, in this strongly 
frustrated approach Tc > Tf (Tc and Tf are the Curie and freezing tempera- 
tures, respectively) in disagreement with experimental results for CeNii-xCux- 
Therefore, the motivation for the present work is to understand better the role of 
disorder as source of frustration in a Kondo lattice model and, for that reason 
to possibly address the experimental findings of CeNii-xCux phase diagram. 

In Fig. 1, the results for a = 0.02 are presented. For high T/J and small 
Jk/J, the numerical solutions display a paramagnetic (PARA) behavior with 
q — 0, m — and A = 0. The solutions remain the same in this small Jk/J 
region until T w 0.61J, where q starts to be continuously non-zero indicating 
a second order phase transition to a spin glass phase. The behavior of the 
order parameters as a function of the temperature is shown in Fig. 2. In 
particular, we can see that at the same temperature where q =/= 0, there is also 
a cusp in X = x/ = (q — q) (x is static susceptibility). If the temperature is 
lowered, the results remain yielding a spin glass solution until T « 0.38 J. From 
that point, the parameters m and q become simultaneously different from zero 
similar to the ferromagnetic solution already found in Ref. 0. However, in the 
present case there is an abrupt change in their behaviour indicating a first order 
transition. In fact, for 0.31J < 0.38J, we have found metastable solutions 
with TO 7^ and q ^ which corresponds to the Mattis states. The emergence 
of these metastable Mattis states at temperature Ti can be seen clearly just 
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Figure 4: The phase diagram with T/J versus Jk/J for a — 0.04 show- 
ing the phases SG (spin glass), FERRO (ferromagnetism) and KONDO (the 
Kondo state). In the region SG+Ferro, there is the coexistence between SG and 
FERRO. The dotted hne means the "pure" Kondo temperature. 



by foUowing the free energy (see Fig. 3). In Fig. 1, the corresponding region 
has been named as SG + Ferro. Finally, at T2, the spin glass solutions become 
unstable thermodynamically while the solution with to 7^ and q 7^ (the 
Mattis states) becomes stable. 

When Jk/J is enhanced (see Fig. 1), it is found a line transition Jk^{T/J) 
for the Kondo state. This kind of solution had already been found in Refs. 
[3 01 E] . The nature of this line transition is complex. It is second order at high 
temperature changing to first order at low temperature. However, there are 
evidences indicating that this complexity could be nonphysical, in fact, related 
to the approximations made in the present approach |28| . 

If the parameter a is increased (for instance a = 0.04, in Fig. 4), we have, 
basically, a phase diagram displaying the same situation already shown in Fig. 
1. However, the spin glass stability range is increased. Finally, for a = 0.15, 
there is no more Mattis states as solution, the spin glass solution is entirely 
dominating for Jk < Jkc{T) in which Jkc{T) is the phase boundary of the 
Kondo state. In this limit, it is recovered the results obtained in Ref. i4j. 

The experimental results for CeNii-^Cux can be now addressed. The phase 
diagrams obtained with bulk methods and /zSR technique 10 show that 
the chemical substitution of Ni produces in the alloy a complex interplay be- 
tween Kondo effect and magnetism where, for example, the spin glass phase is 
always found at higher temperature than the ferromagnetic one (whatever the 
nature of the spin glass region is). Moreover, it also shows the increase of the 
freezing temperature with the decrease of Cu content in the alloy until x ~ 0.4. 
On the other hand, for x < 0.2 there is a considerable reduction of Ce localized 
magnetic moments due to the Kondo effect. 

If it is assumed that the substitution of Cu by Ni can be related to the 
parameters Jk and a, some important aspects of the experimental scenario can 
be reproduced, for example: 

i) For Jk < Jkc(T) and small a (weak frustration). In this regime the mag- 
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Figure 5: The phase diagram with T/J versus Jk/J showing the phase SG 
(spin glass). The dotted hne means the "pure" Kondo temperature. 



netic solutions are dominant with no trace of a Kondo solution. However, the 
combination of high temperature and the local complex randomness prevents 
any kind of stable alignement of the spins (represented by the order parameter 
m). On the contrary, the solutions for the order parameter show a continous 
direct transitions from paramagnetism to spin glass. Eventually, at lower tem- 
perature, the randomness is not enough to keep avoiding any stable alignement 
of the spins. Therefore, the Mattis states start to appear first as metastable 
solutions until becoming the thermodynamical stable solutions of the problem. 
As consequence, in this weak frustration limit, the correct order for the tran- 
sition temperatures is recovered as compared with the experimental situation 
when 0.8 ^x>^ 0.4, where x is the content of Cu [H]. 

As long as the level of frustration a is increasing, the spin glass component of 
the internal field becomes dominant in a larger range of temperature. Therefore, 
it is obtained an enlargement of the spin glass region (see Figure 1, 4 and 5) 
which precedes the onset of the ferromagnetism. As a consequence, the freezing 
temperature Tf is also increased, which resembles the experimental situation 
when a; ^ 0.4 

ii) When Jk > Jkc{T), the magnetic solutions disappear and the Kondo 
state appears as the unique solution in the problem for any value of a which 
means that there is a complete screening of localized spins due to the Kondo 
effect as in the experimental results for the rich Ni region [5]. 

However, there is still some disagreement related to the results obtained 
from ^Sr spectroscopy ^U] and specific heat which suggests the presence of 
nanoclusters which are frozen at some temperature . For lower temperatures, 
those clusters would percolate yielding a ferromagnetism in a similar process 
proposed to explain manganites (HOI. Our results, instead, in the weak frustrated 
limit, indicate a continuous transition to a spin glass phase and, then, at lower 
temperature, the presence of a first order one with a coexistence region between 
spin glass and ferromagnetism (given by the Mattis states). Nevertheless, we 
believe that the present mean field theory is a clear improvement in the sense 
that it provides an effective mechanism which, at least, gives the correct ordering 
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of the magnetic transition temperatures. 

To conclude, in this work, it has been studied the spin glass-ferromagnetism- 
Kondo phase transitions in a Kondo lattice model with a random coupling Jij 
between the localized spins. The Jij is given as a product of random variables 

(/i — i = 1...N). This choice for Jij introduces a parameter a = 

p/N which allows us to control the degree of frustration. Thereby, the balance 
between the two parameters Jk / J and a controls the emergence of the different 
solutions in the problem. For small a (weak frustration) and Jk / J, we have 
found that there is only the presence of spin glass and ferromagnetic solutions, 
and, particularly that the freezing temperature is higher than the transition 
temperature where ferromagnetic solutions are found in good agreement with 
experiment in Ce(Ni,Cu) alloys. For large Jk/J, there is only a Kondo state, 
whatever the value of a is. The results obtained here are interesting for the 
study of the role of disorder, which seems to be better described by taking into 
account an average of discrete values ^'s, rather than a direct average of the 
intersite J values, as suggested also by new experimental results in CeNii-^Cux 
[TT' alloys. 
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